Simulation of the Spherically Symmetric Stellar Core Collapse, Bounce, and 
Postbounce Evolution of a 13 M Star with Boltzmann Neutrino Transport, and Its 

Implications for the Supernova Mechanism 
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With exact three-flavor Boltzmann neutrino transport, we 
simulate the stellar core collapse, bounce, and postbounce 
evolution of a 13 Mq star in spherical symmetry, the New- 
tonian limit, without invoking convection. In the absence 
of convection, prior spherically symmetric models, which im- 
plemented approximations to Boltzmann transport, failed to 
produce explosions. We are motivated to consider exact trans- 
port to determine if these failures were due to the transport 
approximations made and to answer remaining fundamental 
questions in supernova theory. The model presented here is 
the first in a sequence of models beginning with different pro- 
genitors. In this model, a supernova explosion is not obtained. 
We discuss the ramifications of our results for the supernova 
mechanism. 



I. SUPERNOVA PARADIGM 

Core collapse supernovae are among the most impor- 
tant phenomena in astrophysics because of their ener- 
getics and nucleosynthesis. Beginning with the first nu- 
merical simulations conducted by Colgate and White , 
three decades of supernova modeling have established a 
basic supernova paradigm. The supernova shock wave — 
formed when the iron core of a massive star collapses 
gravitationally and rebounds as the core matter exceeds 
nuclear densities — stalls in the iron core as a result of 
enervating losses to nuclear dissociation and neutrinos. 
The failure of this "prompt" supernova mechanism sets 
the stage for a "delayed" mechanism, whereby the shock 
is reenergized by the intense neutrino flux emerging from 
the neutrinospheres carrying off the binding energy of the 
proto- neutron star • The heating is mediated primar- 
ily by the absorption of electron neutrinos and antineu- 
trinos on the dissociation-liberated nucleons behind the 
shock. This past decade has also seen the emergence of 
multidimensional supernova models, which have investi- 
gated the role convection, rotation, and magnetic fields 
may play in the explosion p|-|lO[] . 

Although a plausible framework is now in place, fun- 
damental questions about the explosion mechanism re- 
main: Is the neutrino heating sufficient, or are multidi- 



mensional effects such as convection and rotation neces- 
sary? Can the basic supernova observable, explosion, be 
reproduced by detailed spherically symmetric models, or 
are multidimensional models required? Without a doubt, 
core collapse supernovae are not spherically symmetric. 
For example, neutron star kicks pd[ | and the polarization 
of supernova emitted light Jl2| cannot arise in spherical 
symmetry. Nonetheless, ascertaining the explosion mech- 
anism and understanding every explosion observable are 
two different goals. To achieve both, simulations in one, 
two, and three dimensions must be coordinated. 

The neutrino energy deposition behind the shock de- 
pends sensitively not only on the neutrino luminosities 
but also on the neutrino spectra and angular distribu- 
tions in the postshock region, necessitating exact multi- 
group (multi-neutrino energy) Boltzmann neutrino trans- 
port. Ten percent variations in any of these quantities 
can make the difference between explosion and failure 
in supernova models P,|l3|]. Past simulations have im- 
plemented increasingly sophisticated approximations to 
Boltzmann transport, the most sophisticated of which is 
multigroup flux- limited diffusion [lj,[l^]. A generic fea- 
ture of this approximation is that it underestimates the 
isotropy of the neutrino angular distributions in the heat- 
ing region and, thus, the heating rate |p^ , p7| . It is im- 
portant to note that, without invoking proto-neutron star 
(e.g., neutron finger) convection, simulations that imple- 
ment multi group flux-limited diffusion do not produce 
explosions QJT^]. Moreover, the existence and vigor of 
proto-neu tron s tar convection is currently a matter of 
debate 
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Wilson |20| | implemented an approximation to Boltz- 
mann neutrino transport by using order-of-magnitude 
parameterizations of the neutrino-matter weak interac- 
tions. His models failed to produce explosions. Core 
collapse simulations that implemented exact Boltzmann 
neutrino transport were completed by Mezzacappa and 
Bruenn J2l],^2| . Following this work, we now present the 
findings of a core collapse, bounce, and postbounce sim- 
ulation. Recognizing the need for more accurate time- 
dependent neutrino transport in supernova models, other 
groups have now developed Boltzmann solvers |ES E5| . 
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II. FOUNDATIONS 

We model the explosion of a 13 M Q star, beginning 
with the precollapse model of Nomoto and Hashimoto 
p6j. The core collapse, bounce, and explosion were simu- 
lated with a new neutrino radiation hydrodynamics code 
for both Newtonian and general relativistic spherically 
symmetric flows: AGILE-BOLTZTRAN. BOLTZTRAN 
is a three-flavor Boltzmann neutrino transport solver 
[ 27, , now extended to fully general relativistic flows 
[29|. In this simulation, it is employed in the 0(v/c) limit 
with 6-point Gaussian quadrature to discretize the neu- 
trino angular distributions and 12 energy groups span- 
ning the range from 5 to 300 MeV to discretize the neu- 
trino spectra. AGILE is a conservative general relativis- 



tic hydrodynamics code [29|,^0|. Its adaptivity enables 
us to resolve and seamless y follow the shock through the 
iron core into the outer stellar layers. 

The equation of state of Lattimer and Swesty |3l]] 
(LS EOS) is employed to calculate the local thermody- 
namic state of the matter in nuclear statistical equilib- 
rium (NSE). For matter initially in the silicon layer, the 
temperatures are insufficient to achieve NSE. In this re- 
gion, the radiation and electron components of the LS 
EOS are used, while an ideal gas of 28 Si is assumed for the 
nuclear component. For typical hydrodynamic timesteps 
(~ .1 millisecond), silicon burning occurs within a single 
timestep for T ~ 5 GK therefore, when a fluid ele- 
ment exceeds a temperature of 5 GK in our simulation, 
the silicon is instantaneously burned, achieving NSE and 
releasing thermal energy equal to the difference in nu- 
clear binding energy between 28 Si and the composition 
determined by the LS EOS. 

We investigated the convergence of the net neutrino 
heating rate as the number of Gaussian quadrature points 
and the number of neutrino energy groups in our Boltz- 
mann simulations were varied, as in Messer et al @. In 
the heating region, the 4- and 6-point rates, the 6- and 
8-point rates, and the 12- and 20-group rates differed 
by at most 5 percent, 3 percent, and 3 percent, respec- 
tively. Moreover, during the course of the important first 
300 ms of our simulation, the maximum variation in the 
total energy is ~ 3 x 10 49 erg, which is a few percent 
of the total energy, and the total lepton number is con- 
served to within a fraction of a percent. Note that the 
numerical uncertainty in the net heating rate (which is 
at most 3 percent in our model) is no larger than the 
uncertainty in the total energy. Therefore, any further 
numerical convergence in the computation of this rate 
would be meaningless. 
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FIG. 1. Radial trajectories of equal mass shells in the iron 
core and silicon layer. We also trace the shock, nuclear burn- 
ing, and dissociation fronts, which carve out three regions in 
the (r, t) plane. A: Silicon. B: Iron produced by infall com- 
pression and heating. C: Free nucleons and alpha particles. 



III. STELLAR CORE COLLAPSE, BOUNCE, AND 
POSTBOUNCE EVOLUTION 

Figure 1 shows the radius-versus-time trajectories of 
equal mass shells (0.01M Q ) in the stellar iron core and 
silicon layer during the first 600 ms of postbounce evolu- 
tion. It also shows the shock and nuclear burning front 
trajectories. At 110 ms after bounce, the shock stalls at 
a radius of 230 km and then recedes for the duration of 
the simulation, and no explosion is launched. The shock 
and burning fronts divide the stellar core and silicon layer 
into three regions: A: Silicon. B: Iron produced by infall 
compression and heating. C: Free nucleons and alpha 
particles. 
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FIG. 2. Three-flavor neutrino luminosities and rms ener- 
gies at a radius of 500 km in the core as a function of time 
during our simulation. 
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In Figure 2, we plot the neutrino luminosities and rms 
energies at 500 km in the stellar core as a function of 
time. The electron neutrino luminosity decreases from 
its early "burst" value of 3.5 x 10 53 erg/s as we enter the 
postbounce accretion phase. (The early electron neu- 
trino burst occurs as the shock passes the electron neu- 
trinosphere in the core, liberating the trapped neutrinos 
behind it that are produced by electron capture during 
stellar core collapse.) In the accretion phase, the elec- 
tron neutrino luminosity reaches a maximum of 3.6 x 10 52 
erg/s and then decreases slowly as the mass accretion rate 
decreases. The electron antincutrino and muon/tau neu- 
trino and antincutrino luminosities rise after a hot, delep- 
tonized "mantle" forms beneath the shock (this is the 
region above the cold, degenerate, unshocked core, and 
includes the neutrinospheres). In this mantle, electron 
antincutrinos are produced by positron capture and all 
three flavors of neutrinos and antineutrinos are produced 
by electron-positron annihilation. The electron antincu- 
trino luminosity reaches a maximum of 3.3 x 10 52 erg/s 
and exhibits the same subsequent decrease with decreas- 
ing mass accretion rate. The muon/tau neutrino and an- 
tineutrino luminosities on the other hand reach values of 
only 2.0 x 10 52 erg/s, and thereafter decrease with time. 
For all three flavors, the rms energies increase with time, 
owing to infall into an increasingly deep gravitational 
well. Relative to the electron neutrino and antineutrino 
rms energies, the muon/tau neutrino and antineutrino 
rms energies are larger, with values between 20-25 MeV: 
the muon/tau neutrinos and antineutrinos interact only 
via neutral currents and therefore decouple deeper in the 
core at higher densities. The electron neutrino rms ener- 
gies lie in the range between 10-20 MeV. Relative to the 
electron neutrino rms energies, the electron antincutrino 
rms energies are slightly larger because the electron an- 
tineutrinos decouple at slightly higher densities: the core 
material is neutron rich; therefore, electron antincutrino 
absorption on protons is reduced relative to electron neu- 
trino absorption on neutrons. 

In Figure 3, we plot the mass density, entropy per 
baryon (in units of fee), electron fraction, and velocity 
as a function of radius for various time slices in our sim- 
ulation. In the velocity profiles, the initial outward prop- 
agation of the shock is evident, as is the subsequent de- 
crease in radius at later times. Note the increasing infall 
velocities below the shock as it recedes, reaching values 
of -6000 km/s at 600 ms after bounce. Despite the de- 
creasing shock radius and failed explosion, there remains 
a heating region behind the shock, where the entropies 
continue to increase. This is evident in the entropy pro- 
files. Nonetheless, conditions remain such that an explo- 
sion does not occur in this model in the first second. 
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FIG. 3. Core baryon density, entropy per baryon, electron 
fraction, and velocity profiles at select times during our sim- 
ulation. 



IV. OUTLOOK 



We have presented results from the simulation of the 
stellar core collapse, bounce, and the first 600 ms of post- 
bounce evolution of a 13 M Q progenitor. Spherical sym- 
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metry was assumed, 0(v/c) Boltzmann neutrino trans- 
port was implemented, and gravity was Newtonian. No 
explosion was obtained. In light of our implementation 
of Boltzmann transport, if we do not obtain explosions in 
models initiated from other progenitors (see also Rampp 
and Janka |25|| ), it would indicate that improvements in 
our initial conditions (precollapse models) and/or input 
physics are needed, and/or that the inclusion of mul- 
tidimensional effects such as convection, rotation, and 
magnetic fields are required ingredients in the recipe for 
explosion. In the past, it was not clear whether failures 
to produce explosions in spherically symmetric models 
were the result of transport approximations or the ne- 
glect of an important physical effect. We will report on 
the general relativistic case ^9[[53) and other models in 
subsequent papers. 

Potential improvements in our initial conditions and 
input physics include: improvements in precollapse mod- 
els [p4] p7| ; the use of ensembles of nuclei in the stellar 
core rather than a single representative nucleus; deter- 
mining the electron capture rates on this ensemble with 
detailed shell model computations p^j; the inclusion of 
nucleon correlations in the high-density neutrino opaci- 
ties p9|p0| ; and the inclusion of new neutrino emissivities 
in dense matter [[ll]. These improvements all have the 
potential to quantitatively, if not qualitatively, change 
the details of our simulations. Thus, it is important to 
note that our conclusions are drawn considering the ini- 
tial conditions and input physics used. 

ACKNOWLEDGMENTS 

A.M. is supported at the Oak Ridge National Labora- 
tory, managed by UT-Battelle, LLC, for the U.S. Depart- 
ment of Energy under contract DE-AC05-00OR22725. 
M.L. is supported by the National Science Foundation 
under contract AST-9877130 and, formerly, was sup- 
ported by the Swiss National Science Foundation un- 
der contract 2000-53798.98. O.E.B.M. is supported by 
funds from the Joint Institute for Heavy Ion Research 
and a Department of Energy PECASE award. W.R.H. 
is supported by NASA under contract NAG5-8405 and 
by funds from the Joint Institute for Heavy Ion Research. 
F.-K.T. is supported in part by the Swiss National Sci- 
ence Foundation under contract 2000-53798.98 and as a 
Visiting Distinguished Scientist at the Oak Ridge Na- 
tional Laboratory. S.W.B. is supported by the NSF 
under contract 96-18423 and by NASA under contract 
NAG5-3903. Our simulations were carried out on the 
ORNL Physics Division Cray J90 and the National En- 
ergy Research Supercomputer Center Cray J90 and SV-1. 
We thank the referees for their important questions and 
comments. 



[1] S. A. Colgate and R. H. White, Astrophysical Journal 

143, 626 (1966). 
[2] J. R. Wilson, in Numerical Astrophysics, edited by J. M. 

Centrella, J. M. LeBlanc, and R. L. Bowers (Jones and 

Bartlett, Boston, 1985), pp. 422-434. 
[3] H. H. Bethe and J. R. Wilson, Astrophysical Journal 295, 

14 (1985). 

[4] M. Herant et al., Astrophysical Journal 435, 339 (1994). 
[5] A. Burrows, J. Hayes, and B. A. Fryxell, Astrophysical 

Journal 450, 830 (1995). 
[6] H.-T. Janka and E. Miiller, Astronomy and Astrophysics 

306, 167 (1996). 
[7] A. Mezzacappa et al., Astrophysical Journal 493, 848 

(1998). 

[8] A. Mezzacappa et al., Astrophysical Journal 495, 911 

(1998) . 

[9] C. L. Fryer and A. Heger, Astrophysical Journal 541, 
1033 (2000). 

[10] A. M. Khokhlov et al, Astrophysical Journal 524, L107 

(1999) . 

[11] C. L. Fryer, A. Burrows, and W. Benz, Astrophysical 
Journal 496, 333 (1998). 

[12] J. C. Wheeler, in Cosmic Explosions, A Conference Pro- 
ceedings, edited by S. S. Holt (AIP, Melville, 2000). 

[13] A. Burrows and J. Goshy, Astrophysical Journal Letters 
416, L75 (1993). 

[14] S. W. Bruenn, in First Symposium on Nuclear Physics in 
the Universe, edited by M. W. Guidry and M. R. Strayer 
(IOP Publishing, Bristol, 1993), pp. 31-50. 

[15] J. R. Wilson and R. W. Mayle, Physics Reports 227, 97 
(1993). 

[16] H.-T. Janka, Astronomy and Astrophysics 256, 452 
(1992). 

[17] O. E. B. Messer, A. Mezzacappa, S. W. Bruenn, and 
M. W. Guidry, Astrophysical Journal 507, 353 (1998). 

[18] S. W. Bruenn and T. Dineva, Astrophysical Journal Let- 
ters 458, L71 (1996). 

[19] W. Keil, H.-T. Janka, and E. Miiller, Astrophysical Jour- 
nal Letters 473, Llll (1996). 

[20] J. R. Wilson, Astrophysical Journal 163, 209 (1971). 

[21] A. Mezzacappa and S. W. Bruenn, Astrophysical Journal 
405, 637 (1993). 

[22] A. Mezzacappa and S. W. Bruenn, Astrophysical Journal 
410, 740 (1993). 

[23] S. Yamada, H.-T. Janka, and H. Suzuki, Astronomy and 
Astrophysics 344, 533 (1999). 

[24] A. Burrows et al, Astrophysical Journal, 539, 865 

(2000) . 

[25] M. Rampp and H.-T. Janka, Astrophysical Journal 539, 
L33 (2000). 

[26] K. Nomoto and M. Hashimoto, Physics Reports 163, 13 
(1988). 

[27] A. Mezzacappa and S. W. Bruenn, Astrophysical Journal 
405, 669 (1993). 

[28] A. Mezzacappa and O. E. B. Messer, Journal of Compu- 
tational and Applied Mathematics 109, 281 (1998). 

[29] M. Liebendorfer, Ph.D. thesis, University of Basel, Basel, 
Switzerland, 2000. 

[30] M. Liebendorfer and F.-K. Thielemann, in Nineteenth 
Texas Symposium on Relativistic Astrophysics, edited by 



E. Aubourg, T. Montmerle, J. Paul, and P. Peter (Else- 
vier Science B. V., Amsterdam, 2000). 
J. Lattimer and F. D. Swesty, Nuclear Physics A535, 
331 (1991). 

W. R. Hix and F.-K. Thielemann, Astrophysical Journal 
511, 862 (1999). 

M. Liebendorfer et al, Physical Review D, submitted 
(2000). 

G. Bazan and W. D. Arnett, Astrophysical Journal 496, 
316 (1998). 

H. Umeda, K. Nomoto, and T. Nakamura, in The First 
Stars, edited by A. Weiss (Springer, Berlin, 2000). 

A. Heger, N. Langer, and S. Woosley, Astrophysical Jour- 
nal 528, 368 (2000). 

A. Heger , K. Langanke, G. Martinez-Pinedo, and S. 



Woosley, |astro-ph/0007412 



K. Langanke and G. Martinez-Pinedo, Nuclear Physics 
A673, 481 (2000). 

A. Burrows and R. F. Sawyer, Physical Review C58, 554 
(1998). 

S. Reddy, M. Prakash, J. M. Lattimer, and J. A. Pons, 
Physical Review C59, 2888 (1999). 

T. A. Thompson, A. Burrows, and J. E. Horvath, Phys- 
ical Review C62, 035802 (2000). 



5 



